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Abstract 

A version of scattering theory that was developed many years ago to 
treat nuclear scattering processes, has provided a powerful tool to study 
universality in scattering processes involving open quantum systems with 
underlying classically chaotic dynamics. Recently, it has been used to 
make random matrix theory predictions concerning the statistical proper- 
ties of scattering resonances in mesoscopic electron waveguides and elec- 
tromagnetic waveguides. We provide a simple derivation of this scattering 
theory and we compare its predictions to those obtained from an exactly 
solvable scattering model; and we use it to study the scattering of a par- 
ticle wave from a random potential. This method may prove useful in 
distinguishing the effects of chaos from the effects of disorder in real scat- 
tering processes. 



1 Introduction 

Interest in the dynamical properties of open quantum systems at mesoscopic 
and atomic scales has lead to a rebirth of a form of scattering theory which 
was originally developed to deal with very complicated nuclear collision pro- 
cesses. The origin of the scattering theory that we consider here came from the 
recognition that a collision between two nucleons, one of which is very heavy, 
can lead to the creation of an unstable, but very long-lived, compound nucleus 
which eventually decays. In the late 1930's, Kapur and Peierls 0] used this 



1 



fact to formulate a nonperturbative approach to scattering theory in which the 
compound nucleus was viewed as a stable object which was made unstable by 
weak coupling to the continuum. In the late 1940's, Wigner and Eisenbud || 
developed an alternate version of the Kapur-Peierls theory which lead to the 
very widely used R-matrix approach to scattering theory ||. The idea behind 
both these theories is to decompose configuration space into an internal region 
(reaction region) and an external asymptotic scattering region. As we shall see, 
this approach can be made clean and rigorous if there is an abrupt (in con- 
figuration space) transition between the scattering region and the asymptotic 
region. The internal region can be modeled in terms of a complete set of states 
with fixed boundary conditions on the surface of the internal region. The in- 
ternal states can then be coupled to the external asymptotic states. Bloch || 
and Feshbach f| both showed that a consistent theory requires that the cou- 
pling between the internal and external regions must be singular. In the 1960's, 
Weidenmuller developed a phcnomcnological Hamiltonian approach to nuclear 
scattering based on this picture The Hamiltonian for the interior region 
was based on the shell model of the nucleus. The Hamiltonian for the exterior 
scattering region described the asymptotic states. The strength of the coupling 
between the interior and exterior regions was an unknown input parameter. 

The Weidenmuller approach to scattering theory created a framework with 
which to study the manifestations of chaos in the scattering properties of nu- 
clear systems, as well as mesoscopic and atomic systems. It is now well known 
that chaos manifests itself in bounded quantum systems by affecting the sta- 
tistical distribution of spacings between energy levels |7). For open quantum 
systems, in regimes where scattering resonances are not strongly overlapping, 
underlying chaos affects the statistical distribution of resonance spacings and 
resonance widths Q . If the Hamiltonian of the interior region is chosen to be a 
random matrix Hamiltonian, then predictions can be made as to the statistical 
properties of resonance spacings and resonance widths. These predictions can 
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then be compared to nuclear scattering resonances, or to resonances in elec- 
tron waveguides |l3|],and [|11] and electro-magnetic resonators [[Q] and 
[ p~4| whose cavities have been constructed to yield classically chaotic behavior. 
The agreement between random matrix theory predictions and numerical and 
laboratory experiments is quite good. 

The Weidenmuller Hamiltonian is phenomenological and does not have in- 
formation about the strength of the interface coupling. The random matrix 
theories treat the interface coupling as an input parameter. In a real dynami- 
cal system, however, one needs a way to rigorously determine the coupling at 
the interface of the internal and external regions. A very useful way to fix 
the coupling at the interface was provided by Pavlov . The idea of Pavlov 
was that the surface coupling could be fixed by the requirement that the total 
Hamiltonian (including interior and exterior regions) be Hermitian. 

In this paper, we use a simple textbook scattering problem to illustrate and 
clarify many issues concerning this approach to scattering, and we then use it 
to study the effect of disorder on the scattering process. 

In Section (2), we describe the "textbook" scattering problem and obtain 
exact expressions for the reaction function, the scattering function, and the 
Wigner delay times. In Section (3), we develop the Hamiltonian for interior 
and exterior regions of the scattering system. In Section (4) we derive the 
Hermiticity condition. In Section (5), we derive the reaction function. In Section 
(6), we derive the scattering function and locate resonance energies for the case 
of a smooth scattering potential. In Section (7), we use this theory to study 
the scattering of a particle wave from a random scattering potential. Finally, in 
Section (8), we make some concluding remarks. 
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2 Exact Solution 



We will first consider the scattering of a particle of mass, m, due to the potential 
shown in Fig. (1). The quantum particle enters from the right with energy E 
and is reflected back to the left by an infinitely hard wall located at x = 0. A 
barrier of height, Vo, is located < x < a. The Schroedinger equation, which 
describes propagation of a particle wave, \P(x, t), for all times, t, is given by 

dV(x,t) h 2 8 2 ^{x,t) . \ t / \ 

where h is Planck's constant and the potential, V(x) = oo for (— oo < x<0), 
V(x) — Vo for (0 < x<a), and V(x) — for (a < x < oo). Since the potential 
is infinite for — oo < x < 0, no wavefunction can exist in that region of space. 

The Schroedinger equation for energy eigenstates for Region I (0 < 

x < a) in Fig. (1), is 

-£^ + «^> = ^>- <") 

Energy eigenstates in Region I have the form 

= Asin(fc'x) (2.2) 



where k' = y jpr{E — Vo) and A is a normalization constant. The Schroedinger 
equation for energy eigenstates Q^ix), for Region II (a < x < — oo) in Fig. 1, 
is 

Energy eigenstates in Region II have the form 

&j*(x) = B(e' lkx - S(E)e lkx ), (2.4) 

where k — jftE. The first term describes the incoming part of the energy 
eigenstate and the second term describes the outgoing part of that state. The 
constant, B, is a normalization constant and S(E) is the scattering function. 
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When Vq — 0, then S(E) — 1 and the incoming wave is reflected from the wall 
at x = with an overall phase shift of 7r since we can rewrite the minus sign 
as e I7r . Thus, S(E), contains information about the phase shift of the scattered 
wave due to the potential barrier, Vo, in the region, < x < a. 

We can determine how the two parts of the energy eigenstate are connected 
by the condition that the wavefunction and the slope of the wavefunction must 
be continuous at x — a. These two conditions can be combined into a statement 
that the " logarithmic derivative" of the wavefunction must be continuous at the 
interface. Thus we must have 



1 a d$ 



R{E) ®Ua) dx 



$^(a) dx 



(2.5) 



at the interface. The function, R(E), is called the reaction function and it will 
play an important role in the following sections. The S-matrix is then given by 

where the reaction function, 

„,„, tan(fc'a) 

R(E) - fc ^ Q } . (2.7) 

Since the scattering function has unit magnitude we can also write it 

s ^ =e ie(E) =e -2ika e ie'(E )j (2g) 

where 0(E) is the phase shift of the scattered wave relative to the case for which 
Vq = 0. It is straightforward to show that 

tanW = rSiF and tan & = kaR W- (2 - 9) 

The time delay of the scattered particle due to its interaction with the barrier, 
Vo, is given by 

This time delay is called the Wigner delay time [pq|]. 
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For the scattering system we consider here, the phase angle, O'(E) is given 

by 

6'(E) = 2arctan[4tan(/3')l , (2.11) 
p' 

where (3 = ka and [3 1 = k'a. The Wigner delay time is given by 

. d9{E) 



2ma 2 ( 1 



h 2 \(3' 2 + /3 2 tan(/3') 2 





-'-) 


A/3 


0j 



2ma 2 , 

~w (2 - 12) 



In Fig. (2. a) we show the phase angle, 6(E), as a function of energy, and in 
Fig. (2.b), we show the Wigner delay time for parameter values, % = 1, m = 1, 
a = 1, and Vq — 10. The peaks in Fig. (2.b) occur at values of the energy where 
the particle wave resonates with the barrier region. The fact that the Wigner 
delay time goes negative for low energies means that when Vb^O, the particle 
can be reflected earlier by the barrier, Vq, than it would be if Vq = 0. 

In Fig. (3), we show that exact energy eigenstates as a function of position 
for a range of energies. At resonant energies, the energy eigenfunctions have 
enhanced probability in the scattering region. These resonances are also often 
called quasi-bound states because they are associated with complex energy poles 
of the scattering function. In Fig. (4), we plot those complex energy poles of 
S(E) which give rise to the first three resonance peaks in the Wigner delay time. 
It is important to note that for scattering systems with higher space dimension, 
complex energy poles can have a profound effect on the scattering properties. 
For example, in two dimensional electron wave guides they can completely block 
transmission in some channels [jicjl- 



3 Scattering Hamiltonian 

Most of the interesting information regarding scattering events is contained in 
the energy dependence of the scattering phase shifts and in the location of quasi- 
bound state poles. However, we generally cannot find those quantities exactly 
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as we did in the previous section. There are a number of techniques for finding 
them if the scattering potential is weak, but not if the scattering potential is 
strong. In subsequent sections, we describe an approach to scattering theory 
which can deal with strong interactions, provided they are confined to a localized 
region of space. As we mentioned in Sect. (1), this approach to scattering theory 
was originally developed to deal with nuclear scattering processes, but recently 
has provided an important tool for studying universal properties of scattering 
processes induced by underlying chaos. We shall use this alternate approach 
to scattering theory to study the system considered in Sect. (2), and a similar 
system with a random scattering potential. In this way, we can compare its 
predictions to the exact results which were obtained in Sect. (2). 

Let x denote the position operator and let |x) denote its eigenstates, so 
that x\x) = x\x) . The position eigenstates satisfy a completeness relation, 
dx \x){x\ = 1, where 1 is the unit operator, and they are delta normalized, 
J^^dx (x\x') = S(x — x'). Since the potential energy is infinite for x < 0, all 
states will be zero in that region. We can divide the completeness relation into 
three parts, N = J^^dx \x){x\, Q = J^dx \x){x\, and P = J^°dx \x){x\, so that 
N + Q + P = 1. However, the operator N acting on any state gives zero because 
the potential energy is infinite in that region. Therefore, we can remove N from 
the completeness relation without changing our final results. Thus, from now 
on we will write the completeness relation as Q + P = 1. 

The operators, Q and P , are projection operators. They have the property 
that Q = Q 2 , P — P 2 , and QP = PQ = 0. This is easily checked by explicit 
calculation. The operator Q projects any state or operator onto the interval, 
0<x<a, while the operator P projects any state or operator onto the interval, 
a<£<oo. In other words, if the state, has spatial dependence, \E , (x)=(x|\l/), 
over the interval (0 < x < oo), then the state (xlQI^) = ^(x) for (0 < x < a) 
and the state (x\P\^) = for (a < x < oo). 
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Inside the Region I, (0 < x < a) in Fig.(l), we define a Hamiltonian, 

H QQ =Q(^-f + V )Q, (3.1) 

where p is the momentum operator and m is the mass of the particle. The Hamil- 
tonian, Hqq, is Hermitian and therefore it will have a complete, orthonormal 
set of eigenstates which we denote as Q\<j>j). We can write the eigenvalue prob- 
lem in the region, < x < a, as HqqQ^j) = \jQ\(j)j), where Xj is the j th 
energy eigenvalue of Hqq and j — 1, 2, ...oo. Because there is an infinitely hard 
wall at x = 0, the eigenstates cj>j(x)=(x\Q\(j)j) must be zero at x = 0. We have 
some freedom in choosing the boundary condition at x = a and we will do that 
later. The completeness of the states, Q\4>j), allows us to write the completeness 
relation, J2jQ\<f>j)(<l>j\Q = Q- Orthonormality requires that (<f>j\Q\<t>j') = Sj,j' ■ 
The Hamiltonian for the region, a < x < oo, can be written 

H PP =P(^-f)P. (3.2) 
2m 

Its eigenvalues are continuous and have range, (0<-E<oo). The eigenvector of 
Hpp, with eigenvalue, E, will be denoted P\E). The eigenvalue equation then 
reads, H PP P\E) = EP\E). 

Any state, I'l'), can be decomposed into its contributions from the two dis- 
joint regions of configuration space as |W) = Q\^) +P\^>). We can expand Q\^) 
in terms of the complete set of states, Q\(j)j) and we obtain, 

|*)=^Q|^)+P|*), (3.3) 

3 

where [3j = (4>j\Q\^). The function (a;|(5|W) must be equal to the function, 
(a;|-P|\E'), at the interface, x = a. In addition, the slopes of these two functions 
must be equal at the interface. 

We couple the two regions of configuration space, at their interface, with 
the singular operator, V — C8(x — a)p. The coupling constant, C, will be 
determined later. Then 

H QP = QVP = C dxi / dx \xi)S(xi-xo)S(x -a) — (xo\, (3.4) 

JO Ja dX 
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and 



d 



dx\ \xq)5(xq — xi)S(xi — a)- — (xi\. (3.5) 



/>oo r-a 

Hpq = PVQ =C dx Q 

J a JQ 

It is useful to remember that J^dx 5(x~a) = |, J^dx 5(x—x ) = 1 if < x < a, 
and f^dx 6(x — x a ) = if a < x - Note also that 

■-n 2 d 2 



and 



Hqq = dx \x) 
Jo 

f-OO 

Hpp — dx \x) 

J a 



2m dx 2 

h 2 d 2 



f Voj (x\ 
(x\ 



2m dx 2 

The total Hamiltonian of the system can be written 



H — Hqq + H 



pp i- Hqp + H P Q. 



(3.6) 
(3.7) 

(3.8) 



The energy eigenstates, \E), satisfy the eigenvalue equation H\E) = E\E). The 
energy eigenstates can be decomposed into their contributions from the two 
regions of configuration space, so that 



\E)='£ 7j Q\<f> j )+P\E), 



where 7j = (<f>j\Q\E). The eigenvalue equation then takes the form 



(3.9) 








H QQ ■ 


■ Hqp \ 
• Hqp 


/7i<3l<M\ 

l2Q\<t>2) 


= E 


/7i0l<M\ 
72QI02) 


\Hpq 


Hpq . 


■ Hpp ] 


V P\E) J 




V P\E) ) 



(3.10) 



This yields a series of equations 

HqqI^H + Hqp\E) - EQWj)^, 

for j = 1, 2, ... and 

H PP \E) + ^HpqWj)^ - EP\E). 



(3.11) 



(3.12) 



Before we can proceed further, we must find conditions for Hcrmiticity of the 
Hamiltonian, H . 
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4 Hermiticity Condition 



Consider the arbitrary states, |^) and The condition for Hermiticity of 
these states is that 

(3|If|tf)-(tf|Jf|3)*=0. (4.1) 

We can decompose the states |^) and |S) into their contributions to the two 
disjoint configuration space regions and write them in the form, 



3 

where = (<fo|Q|*) and 

|S)=Q|S) +J P|S)=^ aj 0|^)+P|S), 



(4.2) 



(4.3) 



where otj — (4>j\Q\^). In the interior region, we have expanded |*) and |3) 
in terms of the complete set of energy eigenstates, Q\<fij), of the Hamiltonian, 
Hqq. To simplify the notation, let |^-)=/3jQ|0j), \^)=ajQ\4)j) , = |* P ), 

and P|3) = |Sp). The states can be written 



/ l€i> \ 
16) 1 


and |*) = 


/ Ml) \ 


\\Sp)J 




\i*p>/ 



(4.4) 



If we note that 



<3|If|tt> = «£ 1 |,<6| ) ...<3p|) 



(Hqq 

#QQ 



H, 
H, 



QP 
QP 



\ ( IV'i) \ 

IV>2> 
\l*p)/ 



\ Hpq Hpq . . . Hpp ) 
= YI&IHqqI^) +Y1&\Hqp\*p) +Y1( E p\ H pq\^ + (Zp\H P p\*p), 

3 3 3 

(4.5) 

then the condition for Hermiticity is 



(3|#|tt) - (*|ff|H)* = £[&|A QO hk> - (^1^16)1 
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£[&|ffQp|*p> - (Tpj\H Q p\Zp)*] + Yll( E p\HpQ\Tpj) - (*p|#pq&>* 



+ [(Ep\Hpp\H!p) - (^p\H P p\Ep)*} = 0. 



(4.6) 



We can now evaluate Eq. (4.6) term by term. Let (x\ipj}=ipj(x), (x|£j)=£j(x), 
{x\Ep)=E P (x) and (x\^p)=^p(x). Then 



dx 



dx 



C 



dx 



?i 2 / ,d* P 
7^ -p(°) 



dV'i 
dx 



dx 



dx 



*p(o) 



dS, 



dx 



0, 



(4.7) 



is the Hermiticity condition. For our subsequent discussion, it is useful to note 
that for the special case when the eigenstate, <j>j (x) , has zero slope on the 
interface, so = 0, the Hermiticity condition reduces to 



(Z\H\*)-(*\H\Z)* = -J2Uj( 



2m 



Hp (a) 



j 

d* P 



d* ; 



dx 



dx 



dx 



*p(o) 



dS ; 



dx 



0, 



(4.8) 



If we let |*) = (E 1 ) and |S) = \E) in Eq. (4.8), the Hermiticity condition 
reduces to 

(4.9) 



c h 



where 5 , £;,p(x) = (x|P|.E). This Hermiticity condition allows us to determine 
the coupling constant, C. If we remember that (x\Q\E) — X^7j<A;( x )> then 
continuity of the energy eigenfunction at x = a requires that 



m 



(4.10) 



Thus, the Hermiticity condition allows us to determine the strength of the cou- 
pling at the interface. 
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5 The Reaction Function 



We can use the formalism derived in Sections (3) and (4) to derive an expression 
for the reaction function, R(E). Let us return to Eq. (3.11). If we multiply by 
{<f)j\Q, use Eq. (3.4), and the normalization condition, (</>j\Q\<j>j) = 1, we obtain 



4 3 dx 



= 0. (5.1) 



7j 



Similarly, if we multiply Eq. (3.12) by (E\P, we obtain 

-2^J a dx ^ x) ^~ + 

/•OO 

= E dx^ EP (x)^ E . P (x), (5.2) 

J a 

where * Bi p(x)=(x|P|S). 

Before we can proceed further, we must decide on boundary conditions for 
our states, <pj{x). We can use any complete basis set in the region, < x < 
a, provided they have the property (/>j(0) = 0. We are free to choose the 
boundary condition at x = a. We have seen in Eq. (4.8) that the Hermiticity 
condition becomes especially simple for the boundary condition, = 0. 

a 

Therefore, that is the boundary condition that we will use here. The states 
with these boundary conditions are given by <j>j{ x ) = \/f sm (^f)' f° r 3 °dd 

j-2 2-2 

(j = 1,3,...), and the corresponding eigenvalues are Xj = + V . The 

states, <f>j(x) and the eigenvalues, Xj, should not be confused with the exact 
energy eigenvalues and eigenstates in Region I. The states, <fij(x), only serve 
to provide us a complete set of states in Region I which can then be used 
to construct the exact energy eigenstates in Region I. With these boundary 
conditions, Eq. (5.2) takes the form 

Let us now show that this description of the system leads to well known 
results. We can derive the reaction function originally obtained by Wigner and 
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Eisenbud. Let us return to Eq. (5.1) and solve for 7^. We find 

fi 2 1 ,<N>e,p 



2m E — Xj 3 ' dx 



If we combine Eqs. (4.9) and (5.4), we obtain 



E.P 



2m^-^ E — Xj dx 

3 J 



(5.4) 



(5.5) 



From Eq. (5.5), we can express the logarithmic derivative of the outside wave- 
function, ^t-^ttt 2 — 1 evaluated at the interface (and therefore the reaction func- 
tion, R(E)), in terms of the inside states. We find 



R(E)= 



dx 



2ma^ E — A., 
3 



For scattering in higher dimensional space or for systems with internal degrees of 
freedom the reaction function becomes a matrix and is then called the reaction 
matrix. 

In Eq. (2.8), we obtained an exact expression for the reaction function. If 
we equate Eqs. (2.8) and (5.6), we obtain 

tan(fc'a) h 2 ^(a^ja) h 2 ^ sin 2 (^JL) 
( ' k'a 2ma^ E - X, m.n? 2^> F ft a ir»(y-l)» „ ' 

j 3 J=l 8a 2 m ^0 

(5.7) 

for j odd. Let ft 2 = 2 "] / ° . Then we can write Eq. (5.7) in the form 



tan(Vfc 2 — K 2 a) _ 2 ^ 1 /r s \ 

Vfc 2 - K*a = ^ A, P _ K 2 _ ZLii^lil ( ■ } 

.7 = 1 4a 
odd 

From Ref. [0] one can show that this is just the definition of the series 
expansion of the tangent function in terms of its argument. 

The reaction function can also be expressed as a product of matrices. Let 
us consider the first N eigenvalues and eigenvectors of the the Hamiltonian, 
Hqq (we later let N— >oo). We then write the Hamiltonian, Hqq, in the matrix 
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representation in which it is diagonal (in this representation we denote it as 
H^), and obtain 



(5.9) 



We also introduce the vector of cigenstates of Hqq, but evaluated at the inter- 
face, 

/ Ma) \ 
' Ma) ' 







. 


' ° \ 


frN 

H in - 





A 2 . 


. 




V o 


. 


• X N J 



v N 



2ma 



(5.10) 



\(t) N {a)J 

Then the reaction function can be written in the compact form 



R N (E) = v] v 



El 



N 



H, 



A 



V N = 



2 N 

E" 



2ma 



E-Xi 



(5.11) 



where Iat is the NxN dimensional unit matrix and denotes the Hermitian 
adjoint of %. 

6 The Scattering Function 

It is useful to express the scattering function, S(E), in matrix notation, because 
we can then find a very interesting expression for the complex energy poles of 
S(E). In Eq. (2.7), we obtained an expression for the scattering function in 
terms of the reaction function. In terms of the truncated reaction function it is 



S N (E) =e 



where 



Rn(E) = 



ma 



-2ika 



N 

■E- 



1 + ikaR N (E) 
l-ikaR N {E) 



Bin 2 (^) 



(6.1) 



(6.2) 



3 =i E Vo 

In Table (1), we give the positions of the three lowest energy quasi-bound state 
poles found by using the exact expression for R(E) given in Eq. (2.8). These 
values are exact to the number of digits shown. Also in Table (1), we give the 
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Table 1: comparison of pole positions with exact result 



n 


First Pole 


Second Pole 


Third Pole 


250 


17.9812 -5.0071 


45.0729 -17.6046 


93.2048 -34.5167 


350 


17.9771 -4.9999 


45.0197 -17.5295 


92.8700 -34.2099 


500 


17.9742 -4.9945 


44.9806 -17.4736 


92.6336 -33.9798 


1000 


17.9707 -4.9883 


44.9359 -17.4089 


92.3725 -33.7124 


exact 


17.9672 -4.9821 


44.8921 -17.3448 


92.1262 -33.4466 



positions of those poles found by using Eq. (6.1) for four different truncations, 
N = 250, N = 350, N = 500, and N = 1000. The values obtained using Eq. 
(6.1) converge slowly to the correct answer. 

The scattering function can be written in several other forms as well. Let 
us introduce the column vector, 

/ Ma) \ 

wn — VkavN = 



2m 



02(a) 
\cj) N (a)/ 



Then the scattering function becomes 

S N (E) -- 
where 

K N =w^ N 



liha 


'l + iK N ~ 


^— liha 




l-iK N _ 





1 + 



UK 



N 



1-iK 



N 



E1 N - Hi 



N 



W N . 



Using Eq. (6.4), we can write (we suppress the index N) 
K 1 , 1 

n I - 



1 — %K 1 — iw^ 



n=0 

= w\El-H in ) 



El-Hi. 
1 



-w El - H in 
1 



El — Hs, 



-ww 



El - H. 



■w) 



1 



1 



^El - H in - iww^ J El- H, _ 
The scattering function, Sm{E), can now be written in the form 

M 



w, 



S N {E) = l + 2iw] f {El N -H£) 



(6.3) 



(6.4) 
(6.5) 



(6.6) 



Det[El N - Hg - iw N w^ N ] E1 N - H 



w N . 



(6.7) 
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where M is the adjoint of the matrix, SI — Hi n — iww\ and Det denotes its 
determinant. 

The poles of the scattering function arc given by the condition 

T>et[El N - H** n -iw N w ] N ]=0. (6.8) 

Eq. (6.8) gives an N th order polynomial in S whose solutions are the complex 
energies which locate the poles. At first sight, it would appear that the complex 
energy poles of the scattering function are simply given by the eigenvalues of 
the the non-Hermitian matrix, H^^=H-^ + iwjvwjy However, the column 
matrices, wn depend on energy, E, and therefore the eigenvalues, fii(E), of 
H^j: also depend on energy, E. 

In Fig. (5), we locate the zeros of Det [SI — H e ff] — Det [Si — H e ff — iww j: ] 
in the neighborhood of the first three resonance energies. The zeros satisfy the 
equation, 

Det[Sl-ff e// ] = (S- Ml (S))(S- M2 (S))x...x(S- i u J v(S)). (6.9) 

In Fig. (5. a) we plot (S — /xi(S)) versus S for complex values of S in the 
neighborhood of the first pole. We see that (S — /Lti(S)) = at the energy of 
the first pole. In Figs. (5.b) and (5.c), we plot (S - ^i(S)) and (S — fJ-i(E)) at 
the energies of the second and third poles, respectively. Again we see that they 
go to zero at the energies of their respective poles. 

7 Random potential 

The potential used in Section (2), is a smooth step function which allows us 
to solve the scattering problem exactly. It also allows us to obtain exact ex- 
pressions for the basis state energies, Xj and coupling constants, vn- For our 
simple scattering problem, matrix elements of vn simply alternate between two 
constant values. 

We now use these same methods to study scattering from a random step 
potential. The random potential we use is a sequence of 10 tent-like shapes 
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(upright or inverted) on the interval < x < a with a = 100. We again choose 
h = 1 and m = 1. We can express the potential in the form 

v(x) = io+Y^ ( ^ iq 31 ^ ( x - 1Q J + 1Q ) + ^J- 1 ) Q{wj-x)e(x-ioj+io) 

(7.1) 

The values of Vj are chosen at random from a uniform distribution over the 
interval, — 0.5<Vj< + 0.5. Q(x) is the heaviside function. One hundred different 
realizations of this potential are shown in Fig. 6. 

For this larger value of a, there are about 10 to 12 scattering resonance peaks 
in the energy interval 10.5<£'<11. To hnd the eigenvalues and eigenfunctions for 
this random potential we used an implicit finite difference method. When there 
is no randomness (vj = for all j), the Wigner delay time can be calculated 
exactly as we showed in Section 2. Fig. (7) shows the Wigner delay times when 
vj = for all j. In Figure (7) , the exact result is given by the solid line, and 
the approximate result (obtained by using the reaction matrix series expansion, 
Eq. (5.6)) is given by the discrete dots. The agreement is excellent. In Fig. (8), 
we plot the values of ujy which are obtained from the 100 hundred realizations 
of the random potential when a = 100. The peak points (positive and negative) 
on the oscillating solid line give values of vn for the case Vj = for all j. We 
see that the higher eigenmodes are not affected by the random potential. 

In Fig. (9), we show the Wigner delay times for 100 different realizations 
of random potential. The positions of resonance peaks as well as their widths 
change with different potentials. In Fig. (10) we plot a histogram of values of 
the Wigner delay times shown in Figure (9). This is similar to distribution one 
can get from RMT calculations, except that our distribution has a longer tail 
(see |nj |lj|). Finally, in Fig. (11) we show the distribution of Wigner delay 
time resonance widths. We find that the distribution of delay time widths is 
fairly symmetrically distributed around its average value. 

As we can see, this approach allows us to compute the scattering properties 
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of a variety of shapes of scattering potential with great accuracy. The only 
constraint is that the scattering potential must occupy a well defined region of 
space. 
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8 Conclusions 



The configuration space scattering picture, described in Sections (3)-(6) is a 
basis for Weidenmullcr's phcnomcnological Hamiltonian theory of nuclear scat- 
tering, and it has provided a framework to look for universal behavior in the 
scattering properties of open quantum systems whose underlying classical coun- 
terparts are chaotic [§]. There is now a large body of work [J?|] which shows 
that bounded quantum systems, whose underlying classical dynamics is chaotic, 
have energy eigenvalues and energy eigenfunction which have statistical proper- 
ties similar to those of certain types of random Hamiltonian matrices. Similar 
questions are being asked about the scattering properties of open quantum sys- 
tems For open systems, one can look at the statistical properties of the spacings 
and widths of scattering resonances, and of Wigner delay times. 

The theory described in Sections (3)-(6) can be used to compare the results 
of laboratory and numerical experiments to random Hamiltonian matrix pre- 
dictions. For example, for scattering systems in which no magnetic fields are 
present, the Hamiltonian, Hi n , which is formed from the eigenvalues of Hqq, 
is replaced by a Hamiltonian, H™, which is formed from the eigenvalues of 
a Hamiltonian matrix whose matrix elements are random numbers which are 
Gaussian distributed. The elements of the coupling matrix, W, are also chosen 
from a random distribution. The strength of the coupling constant, C, is left 
as a variable parameter in these random matrix theories. The random matrix 
theory predictions are in qualitative agreement with numerical experiments and 
some laboratory experiments. 

However, a number of assumptions underly the random matrix theory predic- 
tions, and there are a number of issues that remain open. For actual scattering 
systems, the coupling constant, C, is fixed by the detailed dynamics using the 
Hermiticity condition, but in random matrix theories it is a variable parameter. 
In random matrix theories, the coupling matrix, w, is chosen from a random dis- 
tribution. How does it actually look for a deterministic system with underlying 
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classical chaos? This has never been studied. 

Random matrix theories also neglect the dependence of H e f / on the energy, 
E. The actual error made in making this assumption needs to be understood. 
For two dimensional electron waveguides, there is an additional problem that 
at the threshold energy where a new channel opens, quasi-bound states inside 
the cavity can extend far down the leads. |^(| It is not clear how well the 
theory presented in Sections (3)- (6) and therefore the random matrix theory, can 
describe scattering in those regions. Never-the-less, this approach to scattering 
theory has allowed contact to be made between chaos theory and scattering 
theory for open quantum systems. As we have just shown, it can also provide 
a powerful tool to study scattering processes in atomic and mesoscopic devices, 
when intrinsic disorder in the medium needs to be included. 
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Figure 1: The scattering potential, V(x) for U = 1, Vq = 10, and a = 1. 



Figure 2: (a) Phase 9 of S = exp(i6) function, (b) Wigner delay time r Both 
plot is as a function of energy for Ti — 1, Vb = 10, and a = 1. Three resonances 
are shown with arrows at the resonance energy values. 

Figure 3: Absolute value of exact energy eigenstates plotted as a function of 
spatial variable, x, and energy, E, for % = 1, Vb = 10, and a = 1. 

Figure 4: Poles of S function in the lower complex plane for fi = 1, Vq = 10, 
and a = 1. The absolute value of S is plotted in logarithmic scale, as a function 
to real and imaginary parts of the energy. 

Figure 5: (a) E — H\(E) versus E = E r + iEi in the neighborhood of the first 
pole; (b) E — ^{E) versus E = E r + iEi in the neighborhood of the second 
pole; (c) E — Hz{E) versus E = E r + iEi in the neighborhood of the third pole. 
All plots are for fi = 1, Vb = 10, and a = 1. 

Figure 6: One hundred realizations of the random potential, Eq. (7.1). 

Figure 7: Wigner delay time for the step potential, V{x) = 10 for < x < 100 
and V(x) — for x > 100. The solid line is the exact result. The dots are 
obtained using the series approximation to the reactions matrix, Eq. (5.6) 

Figure 8: The coupling constants, vjy, for the 100 different realizations of 
random potential shown in Fig. (6). The peak values (positive and negative) of 
the solid line give the values of vn for Vj = and V(x) = 10. 

Figure 9: Wigner delay times for the 100 realizations of the random potential 
shown in Fig. 6. 

Figure 10: The normalized probability distribution, P(t), of scaled Wigner 
delay times, r/(r), where (r) = 125.42 is the average value of r, taken over all 
100 realizations of the random potential. 

Figure 11: The normalized probability distribution, P(T), of scaled Wigner 
delay time half widths, r/(r), where (r) = 0.05 is average half- width taken 
over all 100 realizations of the random potential. 
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